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ABSTRACT 


We  consider  the  dynamo  action  produced  by  convection  of  a partially- 
ionized,  electrically  conducting  gas  in  a magnetic  field.  The  model 
consists  of  two  thin,  Cartesian  unipolar  inductors  connected  in  series 
by  the  magnetic  field.  For  the  case  of  a uniform  magnetic  field  we 
compute  the  total  current  system  generated  by  an  arbitrary  gas  flow; 
for  the  case  of  a non-uniform  field,  we  compute  only  the  field-aligned 
coupling  current.  Application  is  made  to  the  solar  atmosphere. 
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I . INTRODUCTION 


Magnetic  fields  play  a significant  role  in  a wide  variety  of  plane- 
tary, interplanetary,  and  astrophysical  contexts.  One  important  possi- 
bility is  that  plasma  motion  across  field  lines  can  serve  as  a dynamo 
to  convert  mechanical  energy  into  electric  and  magnetic  energy.  Dynamo 
effects  are  conspicuous  in  the  Earth's  magnetosphere  and  ionosphere, 
Jupiter's  magnetosphere,  and  on  the  solar  surface.  Accurate  modelling 
of  these  processes  will  advance  our  understanding  of  solar  system  plasma 
physics  a great  deal. 

Solar  active  regions  are  known  to  be  permeated  by  strong  (~  kG) 

magnetic  fields  in  the  photosphere.  The  large  energy  output  of  flares 

is  generally  attributed  to  the  rapid  conversion  of  magnetic  energy  of  non- 
— ♦ . 

potential  (V  x B r 0)  fields  in  the  corona  into  particle  energy.  In  the 

2 

low  3 = nkT/B  coronal  plasma  above  a well -developed  active  region,  field- 
aligned  currents  can  flow  which  close  in  the  photosphere  and  lower  chromo- 
sphere. The  electromotive  force  for  these  coronal  currents  can  be  under- 
stood as  arising  from  the  convective  motions  of  the  partially  ionized  gas 
in  the  photosphere  at  a pressure  p ~ 0.1  atm.  If  we  accept  a value  of 
7]  ~ 10  ohm  cm  for  the  intrinsic  resistivity  of  the  photosphere,  then  the 

gas  has  electrical  conducting  properties  equivalent  to  a semi-conductor 
with  a conductivity  somewhat  better  than  sea  water.  Photospheric  "neutral 
winds"  [Sen  and  White,  1972]  can  blow  across  the  magnetic  field  and  impart 
stresses  on  it  causing  it  to  be  twisted  or  sheared.  The  flare  is  then  a 
rapid  release  of  this  free  energy. 

The  simplest  approach  to  the  dynamo  problem  is  to  assume  infinite 
conductivity  in  the  photosphere  for  which  the  field  lines  are  frozen  into 
the  gas.  For  time  scales  less  than  the  resistive  time  constant,  this  is 
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probably  an  accurate  picture.  Finite  resistivity  is  generally  incor- 
porated by  solving  a convection/diffusion  equation  for  the  magnetic  field. 

However,  there  are  advantages  to  an  alternate  perspective  of  the 
problem  given  by  a source-theoretic  analysis.  Here,  one  solves  directly 
for  the  currents  using,  say,  Ohm's  law  and  self -consistently  adjusts  the 
magnetic  field.  This  approach  has  the  advantage  of  enabling  one  to  assess 
and  incorporate  anomalous  plasma  and  kinetic  effects  (e.g.,  current -driven 
plasma  instabilities)  and  also  to  handle  tricky  boundary  conditions  and 
constraints  which  can  arise  when  the  total  system  consists  of  quite 
different  behaving  plasmas.  Alfven  and  Falthammer  (1963)  have  described 
three  classes  of  plasmas  with  radically  different  properties  and  when 
combinations  of  these  occur  in  a system,  source  theory  is  well  suited  for 
the  analysis  on  a local  as  well  as  global  basis. 

Sen  and  White  (1972)  adopted  the  picture  of  the  photosphere  as  an 
MHD  dynamo  and  considered  the  currents  driven  by  radial  gas  flow  from  the 
center  of  a symmetric  sunspot.  They  only  considered  those  currents  which 
flow  in  the  photosphere  and  their  analysis  could  not  describe  any  twisting 
or  shearing  of  the  field.  Heyvaerts  (197*0 > however,  provided  the  key 
idea  that  the  magnetic  field  could  act  as  a shunt  and  allow  a large  current 
flow  out  of  the  photosphere  and  into  the  corona.  This  aspect  character- 
izes the  photosphere  as  a full  MHD  generator  [Brogan,  1968]  capable  of 
also  storing  magnetic  energy  through  a coronal  current  system. 

This  picture  is  attractive  for  flares  because  if  photospherlc  faculae 
and  chromospheric  plages  are  interpreted  as  resulting  from  Joule  dissipa- 
tion of  dynamo  currents  [Sen  and  White,  1972],  then  in  view  of  Svestka's 
(1976)  description  "A  flare  appears  either  as  a brightening  of  a part  of 
existing  plage  (the  most  common  case)  or  as  the  formation  of  new  bright 


areas  in  places  where  the  plage  did  not  exist  before",  the  plage  is  a 
signature  of  dynamo  action  (either  as  a generator  or  a load)  and  its 
coronal  magnetic  field  can  be  highly  non-potential  and  susceptible  to  a 
flare  instability. 

The  purpose  of  this  paper  is  to  synthesize  the  work  of  Sen  and 
White  (1972)  and  Heyvaerts  (197*0  in  a model  which  allows  more  general 
gas  flows.  The  very  low  electrical  resistance  of  the  corona  compared  to 
the  photosphere  suggests  a very  simple  geometrical  model  of  the  system 
consisting  of  2 + 1 dimensions  can  be  used.  We  consider  the  dynamo  action 
in  a relatively  thin  (100-1000  km)  horizontal  layer  on  the  top  of  the 
photosphere  and  allow  currents  to  pass  with  little  resistance  into  the 
corona  and  dissipate  in  a photospheric  load  at  the  conjugate  magnetic 
footpoints.  The  model  can  be  viewed  as  simply  two  resistive  unipolar  in- 
ductors connected  in  series. 

The  most  complete  analysis  of  this  problem  would  require  that  we 
self-consistently  compute  the  generated  magnetic  field.  But  for  simpli- 
city we  shall  assume  1)  steady-state  conditions  and  2)  the  magnetic  Reynolds 

number  R < 1 . The  second  condition  puts  us  in  the  regime  of  conventional 
m 

MHD  (electric)  generators  [Brogan,  I968] . However,  the  analysis  for 

R » 1 is  inherently  non-linear  and  our  simplified  linear  treatment 
m 

should  contain  useful  information  even  for  R **<  1 . It  may  well  be 

m 

that  an  MHD  or  tearing  instability  prevents  a very  large  deviation  from 
the  potential  magnetic  field  configuration. 

In  Section  II  we  describe  the  model  of  two  coupled  parallel-plane 
unipolar  inductors.  In  Section  III  we  compute  the  Hall,  Pedersen,  and 
field-aligned  currents  for  an  arbitrary  two-dimensional  neutral  wind 
velocity  flow  for  the  case  of  uniform  B . In  Section  IV  we  consider  an 
inhomogeneous  B and  compute  only  the  field-aligned  current.  In  Section  V 
we  present  concluding  remarks. 
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II.  MODEL 


Heyvaerts'  (197*0  shunt  model  of  the  photosphere  is  shown  schemati- 
cally in  Fig.  1(a)  for  a simple  bipolar  magnetic  field  configuration. 

The  dynamo  slabs  are  linked  by  the  magnetic  field  and  because  of  the  low 
electrical  resistivity  of  the  fully  ionized,  low  p plasma  in  the  corona 
the  field  lines  are  near  iso-potential  contours  of  the  electric  field. 

We  have  modified  his  model  slightly  and  will  consider  the  system  sketched 
in  Fig.  1(b).  This  configuration  simplifies  the  analysis  somewhat;  it 
also  lends  more  generality  to  the  problem  both  in  execution  and  applica- 
tion to  other  space  and  astrophysical  circumstances.  We  are  only 
interested  in  the  current  systems  possible  and  the  generated  magnetic 
fields  may  differ  for  1(a)  and  1(b)  because  of  geometrical  considerations. 
We  have  also  neglected  effects  due  to  horizontal  components  of  the 
magnetic  field  in  the  slabs  for  the  sake  of  simplicity. 

We  assume  a neutral  wind  flows  in  each  slab  of  effective  width 
-L&  b for  which  the  two  dimensional  (anisotropic)  electrical  conductivity 
tensor  o perpendicular  to  the  magnetic  field  is  known.  The  thickness 
of  the  slab  is  treated  as  small  so  that  the  potential  drop  along  the 
magnetic  field  is  small  everywhere.  The  magnetic  field  lines  are  equi- 
potentials  of  the  electric  field  and  conduct  current  between  each  slab 
with  no  resistance;  in  the  interface  region,  the  magnetic  field  lines 
are  perfect  insulators  in  the  perpendicular  direction. 
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III.  HOMOGENEOUS  MAGNETIC  FIELD 


A.  CURRENTS 

If  the  magnetic  field  is  uniform  in  the  dynamo  region,  the 
analysis  yields  differential  equations  with  constant  coefficients  and  a 
complete  solution  for  the  currents  may  be  obtained,  lor  two  dimensional 
cartesian  coordinates  we  write  the  conductivity  tensor  in  a slab  as 


v-v'ik  “He  i) 


which  defines  the  matrices  v and  e . The  scalar  quantities  and 

CTg  are  the  Pedersen  and  Hall  conductivities  respectively.  If  we  can 
neglect  the  effects  of  electron-ion  collisions,  the  conductivities  may 
be  written  as 


CT1  = n0e 


m(v2  M(vT  +q2J 

L v en  ce  v in  ci'- 


(1  + a) 


CT2  = n0e 


m(v2  +n2  ) m(v?  +cP  ) 

_ en  ce'  v in  cl'. 


Here  we  have  assumed  a simple  gas  where  the  charge  carriers  are  electrons 

of  mass  m and  singly  charged  ions  of  mass  M with  densities  n n.  = n^ 

e i 0 

The  approximate  forms  of  (2)  and  (3)  result  for  the  dynamo  inequalities 


[ Sen  and  White,  1972] 


v /n  « 1 , v.  /0  . » 1 

en  ce  * in  ci 


which  can  occur  in  the  photosphere  in  the  presence  of  kG  magnetic  fields. 

VenV 

When  a = zr — pr — < 1 the  ion  slip  regime  of  conventional  MHD  generators 


“ce  ci 

occurs  [Swift-Hook,  1969] • We  also  may  use 


(5) 


and 


= C1  1 + (a2/<Jl)2 


(6) 


for  the  intrinsic  (field-aligned)  and  Cowling  conductivities  in  the 
photosphere.  A thorough  discussion  of  electrical  conductivity  may  be 
found  in  the  reviews  of  Chapman  (19^6)  and  Braginskii  (1965)*  An  excel- 
lent review  of  MHD  generator  physics  can  be  found  in  Swift-Hook  (I965) . 

Referring  to  Fig.  1(b),  if  the  field-aligned  current  density  falls 
off  rapidly  inside  the  slab  and  in  the  vicinity  of  z = 0 , we  may 
approximate  J^(x, y, z)  =-  J^(x, y)H(z)  where  H(z)  is  the  Heaviside  step 
function,  then  integration  of  the  charge  conservation  equation  yields 


so  that 


3J  Sj  \ 

—L  + _L 

3y  Sz  / 


= 0 


(7) 


where 


and 


is  the  surface  current  density  in  the  slab  across  the  magnetic  field. 
All  spatial  dependence  of  the  variables  now  occurs  in  the  (x, y)  plane 
and  we  shall  make  use  of  the  (scalar)  two  dimensional  curl  operation 
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defined  by 


dv  8v 

V x V * V • eV  = 3/  - ^ 


For  the  case  of  a uniform  magnetic  field,  slab  a maps  identically  onto 
b and  we  have 


Jz(x,y)  = - V • Ka(x,y)  = V • K^y)  . 

If  for  each  slab  we  define  the  height-integrated  conductivity  as 


(10) 


/ 


^2  = /dza~CT't, 


(11) 


where  •{,  measures  the  effective  thickness  of  each  slab,  then  Ohm's  law 
may  be  written  in  matrix  form  as 


. =V  [e  + — eV  J . 

a'b  ^a,bL  C ^bJ 


(12) 


E(x, y)  is  the  electric  field  which  in  steady  state  satisfies 

V X E = 0 j 

V • E = 4jtp  ) 


(13) 


and  the  second  term  in  the  brackets  is  the  electromotive  field  resulting 

from  a neutral  wind  velocity  field  V (x, y)  in  each  slab. 

a,  D 

We  now  add  the  expressions  in  (12)  for  a and  b and  using  (10) 
and  (13)  we  obtain 

V*E  = V*  |(s  V + s V_  ) + s (sn  V + s.  V )]  , (14) 

sla+s^b  |_  2a  a 2b  b la  a lb  b J ’ ' 

where  we  have  defined  analogous  to  (1) 


J2  =*  sifc  - s2e 


(15) 
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and  have  used  the  relation  ee  = - t . 

It  is  convenient,  at  this  point,  to  introduce  the  potential  functions, 


ijr  and  x , such  that 


V = V,  + V = - V*  + eVX  . ■ 
t X 


V and  V are  the  irrotational  and  solenoidal  components  of  the  neutral 

i|)  X 


velocity  field  satisfying 


V2*  = 


f = - V • V 


V^X  = - 


V X v /'  . 


Defining  the  electric  potential  to  be  E = - V0  , we  obtain  from  (14) 

**  ■ 77%  * *1/*. + * -nA,]  • (is) 


The  solution  to  (18)  is  given  as 


0 = — — I S \J|  + S-  ' 

sla+slb  L 2*  * la 


.X.  + s2b*b  * slb*b]  + A • 


A is  the  homogeneous  solution  satisfying  V^A  = 0 and  represents  the 

effects  of  any  externally  imposed  boundary  conditions  in  the  (x, y) 

plane.  However,  for  the  geometry  of  Fig.  1(b)  it  is  assumed  that  all 

sources  of  the  electric  field  have  been  accounted  for  by  specification 

of  V . . We  may.  therefore,  take  A = 0 and  the  electric  field  is 

a,  b 

given  as 

E = - - / -C—  ("s  V\|r  + sn  Vx  + s Vijr  + s ^X  1 . (2C 

si  2a  Y a la  a 2b  yb  lb  b v 

la  lb  L J 

Inserting  this  expression  into  (12)  we  have  finally 
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r 


1 


K 

a 


B/c 


s,  + S... 
la  lb 


s Vtlf 
2b  vb 


s-.Vi|r.  + “ sni(s^  V*  .+ 

2b  Yb  lb  b'  lbv  2a  ’a 

p p 

+ S-ik7*^)  + (Si  +So  +Sn  St  w 

lb  b la  2a  la  lb'  fa 


6l.V] 

(21)  | 
- slbs2«ma]  ) 


B/c 


(. 


b s,  + s_ 
la  lb 


|I{sla(s2b'VslbWb)  - WV9!,7*/] 


(22) 

I 

I 


J = 


-efs  (s  7*  +s  VX  ) + (s2  + s2  +s,  s ) V\|r  - s s VX  ] ! 

|_  2bx  2a  fa  la  a'  '’lb  2b  la  lb'  ’b  la  2b  bj  ( 

— — j s (s  V2*  +s  V^x  ) - s ( s„.  V2#  + s.. . V^x.)  I (23) 

3la+s^b  J lbv  2a  T a la  a ' lav  2b  yb  lb  b' | v ' 


We  remark  that  all  ambipolar  effects  are  automatically  taken  account 
of  in  the  solution  for  the  polarization  electric  field  (20).  The 
representations  in  (21)  and  (22)  show  explicitly  the  irrotational  and 
solenoidal  components  of  K in  each  slab,  the  former  giving  rise  to  the 
field-aligned  current  (23)  linking  both  slabs  either  of  which  can  be  a 
"generator/motor"  or  just  simply  a "load". 

I f we  rewrite  (12)  for  each  slab  (suppressing  subscripts)  as 


E + — V X B 
c 


E"1-  K 


(24) 


and  make  use  of  the  fact  that 

K-r"1  , 

«.  ‘3 

then  scalar  multiplication  of  (24)  with  K yields 


(K  X B) 


E . 


(25) 


(26) 
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If  we  multiply  (26)  by  an  infinitesimal  element  of  area  dA  , then 


the  first  term  is  the  work  done  by  the  height-integrated  volume  force 


— * — • — * 

— K X B at  velocity  V and,  if  negative,  represents  the  local  power 


generated  by  the  slab  as  a seat  of  electromotive  force.  The  second  term 
represents  ohmic  losses  and  the  right-hand  side  gives  the  electrical 
power  into  dA  . A negative  value  of  K * E represents  electrical  power 

— » "4  — * 

out  of  dA  and  when  this  occurs  a negative  value  of  V • (K  X b)  is 
required;  the  slab  is  then  locally  a power  generator  [Swift-Hook,  1965] . 

“4  -4  -4  . 

We  note  that  in  the  most  general  case  the  quantity  ( Ka  + 1C  )*E  r 0 
at  a point  (x, y).  This  seems  to  violate  conservation  of  electrical 
energy.  However,  it  is  easy  to  show  that  if  we  integrate  over  the  entire 
area  of  the  system,  J dA  (K&  + K^’E  = 0 . This  enforces  conservation 
of  total  electrical  energy  and  indicates  that  different  points  in  a slab 

— * Q “ ♦ — « 

communicate  power  through  the  Poynting  flux  S =.  ^ E X B . 

B.  DECOMPOSITION  OF  THE  VELOCITY  FIELD 

The  solutions  (20)-(23)  have  been  readily  obtained  because  of 
the  possibility  of  the  decomposition  in  (16).  The  Fourier  expansion  of 

V^G  = - 4jt6(x-x/)6(y-y/) 


(27) 


leads  to  the  Green's  function  in  two  dimensions  (k  = k^  ) 


H(k)  - % ■ 


(28) 


Thus,  for  a given  neutral  velocity  field  in  each  slab,  the  solution  of 

(17)  is 


(*)  f dK  1 ‘W-,  -rV<  (v'*vl 


(29) 


and  the  formal  solution  of  the  problem  is  complete. 
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C.  SIMPLE  CAS  FLOWS 


Consideration  of  some  elementary  gas  flows  in  special  circum- 
stances will  illuminate  much  of  the  physics  contained  in  (20)-(23).  We 
take  Vb  = 0 and  let  = V for  the  following  cases  of  interest: 

(i)  OPEN  TERMINAL  (£  = q) 


Current  flow  is  restricted  to  slab 
The  electric  field  (20)  reduces  to 


a and 


Jz  = \ = 0 


E = - 


2a 

Sla 


* y 
c \|| 


1 -*  — • 
- vv  x B 
c X 


(30) 


and  the  surface  current  (21)  is 


Ka  = c S3a  Vt  * B ' 


(31) 


The  magnetic  force  density  on  the  gas  is  given  as 


1 — * — * 
- K X B 
c 


I2  S3a  V* 


(32) 


while  the  ohmic  losses  and  electrical  power  density  reduce  to 


3a 


V2 

* 


(33) 


K 


E 


2 

c 


3o  V, 

3a  i|r 


(3*0 


The  quantities  in  (33)  and  (3*0  will,  in  general,  differ  when  dA  is  a 
seat  of  EMF.  We  could  "redefine"  the  electric  field  E E to  absorb 
the  V X b term  in  (12) ; then,  of  course,  — K = K • E . However,  we 

S3 

prefer  to  distinguish  between  the  electromotive  field  and  the  electric 
field  satisfying  (13)  which,  in  principle,  can  unambiguously  be  deter- 
mined by  hooking  up  a potentiometer  across  the  slab's  magnetic  field  terminals. 


12 


We  note  that  only  the  irrotational  component  of  V experiences  a 
braking  force  in  (32).  The  solenoidal  part  of  V flows  freely  due  to 
the  absence  of  any  electrical  load  in  slab  b.  We  also  note  that  the 
surface  current  in  (31),  which  is  manifestly  solenoidal,  is  driven  by 
the  irrotational  component  of  V by  means  of  the  "enhanced"  or  Cowling 
conductivity  s^  . Clearly,  a polarization  electric  field  (30)  has 
built  up  which  enhances  the  current  flow  above  that  expected  from  the 
Pedersen  conductivity  alone.  If,  also  2^  = , the  polarization  fijgld 

1 sp  — — 

will  drive,  in  addition  to  other  currents,  a strong  component  = — — V, X B 
’ ' 2c  s^  f 

in  slab  b.  By  analogy  to  the  effect  observed  in  the  earth's  equatorial 

ionosphere  [ Cunnold,  I978]  we  may  consider  the  current  in  (31)  driven 

by  V through  the  Cowling  conductivity  as  a two-dimensional  electrojet. 

* 

Although  the  effects  are  similar,  in  fact  the  earth's  electrojet  (s)  is 
much  more  complicated  owing  to  the  complexities  of  the  three-dimensional 
field/ionosphere  configuration. 


(ii)  SHORTED  TERMINAL  ( slb  - 00  , 


2b 
5 lb 


0) 


Of  course  E = 0 and  the  surface  currents  are 


-»  B “♦  1 — 

K = - s.  V+-s,  V X B 

a c 2a  c la 

-»  b 1 -»  -» 

1C  = S_  V, S,  V X B 

b c2a\|i  clax 


(35) 


It  is  interesting  here  to  compare  the  body  forces 


Fa  = \ S2a  V B ' ^ ‘l«  \ - % *la  \ + \ s2a  \ X 9 

c c 


(36) 


F = - ~ s V,  X B 
b 2 2a  \|r 


+ 2 Sla  VX 

c 
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For  an  irrotational  flow,  the  magnetic  field  opposes  the  gas  flow  and 
slab  b exerts  a force  perpendicular  to  the  gas  flow;  a solenoidal  flow 
is  braked  by  slab  b and  the  magnetic  field  exerts  a force  perpendicular 
to  the  gas  flow. 

( iii)  SOLENOIDAL  FLOW  (v  * V = 0) 

The  currents  are  given  as 


and 


K = — lb 


. - , Bs  V + s V 

« c Sla+Slb  L 2a  X la  X 


x ij 


K , - 1 _!a_  f 

b c sn  +sn. 

la  lb  L 


Bs„  V + s.  V x B 
2b  X lb  X 


(37) 


J 

Z 


B SlaSlb 
C Sla+Slb 


V2 


The  electric  field  is 


(38) 


i.-i 


la 


V X B 


(39) 


C Sla+Slb  X 

and  becomes  E — V.,  X B only  when  either  s, . 0 or  s 

c X lb 

Heyvaerts'  (197*0  model  falls  into  this  class  of  photospheric  flows  which 
are  good  candidates  for  shearing,  say,  a bipolar  magnetic  arcade. 


la 


(iv)  IRROTATIONAL  FLOW  (V  X V = 0) 

Sen  and  White’s  (1972)  model  falls  into  this  class  and, 
in  light  of  recent  solar  observations  of  sub-arcsecond,  kilogauss 
magnetic  flux  elements,  we  reexamine  their  model.  The  currents  are 
given  as 
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K = — |bs  s V + (si  +si  +s.  s 

a c s.,  + s. . I ootkj.  vio  o» 

lfl  lh  *- 


2 2 

|_“”2a“lb'i(r  v^la  S2a  "la"lb' 


>\ x ST 


(40) 


2a 


c s.  + S... 
la  lb 


and 


Bs  V + s V x B 
lb  f 2b  ^ 


j = * S-^-lb-  V2, 

z c s,  +r 


la  "lb 


(*1) 


The  corresponding  electric  field  is 


- B S2a  - 
E V, 

c si.”ib  * 


m 


Irrotational,  radial  gas  flow  may  be  described  by  the  expression 

V = v(p)[h(p-Pi)  - h( p -Pg)J p , 

where  p^  < p^  defines  the  inner  and  outer  boundaries  of  a flow  annulus. 

Sen  and  White  (1972)  considered  an  outflow,  but  we  shall  instead  consider 
an  inflow,  V(p)  < 0 . 

If  we  take  E = 0 , the  above  expressions  reduce  to  the  results  of 
b 

Sen  and  White  (1972)  for  the  ion  slip  regime  they  considered.  In  this 
situation  the  electrons  are  strongly  magnetized  and  the  ions  unmagnetized. 
Physically,  the  gas  carries  ions  inwards  until  a polarization  field  builds 
up  to  cancel  any  radial  electric  current,  = 0 . A small  azimuthal 

— * *4 

ion  current  flows  in  the  V X B direction.  However,  in  response  to 

■4  «4 

the  polarization  electric  field,  electrons  E x B drift  and  enhance 
the  azimuthal  current  considerably.  This  intense  azimuthal  current  system 
(primarily  of  electrons)  was  termed  a Hall  current  by  Sen  and  White  in 
reference  to  the  direction  of  the  electric  field;  however,  in  the  language 

— • -4 

of  MHD  generators,  a current  flow  in  the  V x B direction  is  conventionally 


1 
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a Faraday  or  Pedersen  current  [Swift-Hook,  I965 ; Brogan,  I968]  in 
reference  to  the  direction  of  the  electromotive  field.  It  is  clear  that 
either  convention  is  acceptable  so  long  as  one  makes  explicit  his  refer- 
ence to  either  a driving  electric  or  electromotive  field. 

Since  inward  gas  flow  enhances  a "seed"  magnetic  field,  the  question 
is  can  a self-consistent  magnetic  dynamo  be  set  up.  If  the  electrojet 
current  is  confined  to  a thin  sheet,  it  is  straightforward  to  show  that 
in  steady-state,  the  velocity  required  to  maintain  the  magnetic  field 
must  satisfy  the  order  of  magnitude  estimate  (esu) 


This  result  is  in  agreement  with  Cowling's  (1957)  estimate  (aside  from 
a factor  of  two) ; however,  we  emphasize  that  the  conductivity  to  be 
used  is  the  Cowling  conductivity  (6)  and  the  scale  length  is  the  thick- 
ness of  the  current  sheet  not  the  perpendicular  scale  length,  r . For 
a current  sheet  where  -t  « r the  magnetic  field  does  not  depend  sensi- 
tively on  the  vertical  distribution  of  current  inside  the  sheet  so  long 
as  the  surface  current  density  is  held  constant,  K “ -L-V  . Thus,  the 
relevant  combination  must  be  as  given  in  (4j).  If  we  take  values  of 

~ 1011  esu  and  -L  ~ 100  km.  Equation  (43)  yields  V = 1.5  m/sec.  This 
value  is  10^  times  larger  than  Cowling’s  (1957)  estimate  since  his  conducti- 
vity was  a coronal  value  and  he  used  a very  large  perpendicular  scale 
length. 

The  very  short  lifetime  (^  hrs)  of  sub-arcsecond  magnetic  flux 
elements  (fluxules)  [Harvey,  1977 1 Stenflo,  1976;  Parker,  I978,  and  his 
references]  strongly  suggests  that  they  are  a solar  surface  phenomenon 
and  are  generated  in  the  photosphere  (or  not  far  from  it).  Parker's  (I978) 
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analysis  demonstrated  that  a superadlabatic  effect  occurring  beneath  the 
fluxule  could  lower  the  temperature  of  the  core  and  thereby  reduce  the 
core  gas  pressure.  An  effect  of  this  type  can  produce  an  MHD  magnetic 
generator  where  the  core  is  a low  pressure,  MHD  sump.  Observations  clearly 
show  gas  downdrafts  into  the  core  and  we  may  expect  also  a smaller  hori- 
zontal velocity  component  across  the  field  into  the  core. 


Parker's  analysis  also  raises  an  important  point  concerning  the 

electric  conductivity.  The  fact  that  ion-neutral  collisions  are  important 

in  the  quiet  photosphere  has  already  been  appreciated  by  Kopecky  and 

Kuklin  (1969)  and  by  Sen  and  White  (1972).  In  this  regime  the  conductivity 

is  very  sensitive  to  gas  temperature  relating  directly  to  the  degree  of 

ionization.  A cool  core  may  even  lower  the  conductivity  by  as  much  as 

two  orders  of  magnitude  from  the  value  quoted  earlier  and  place  the  effec- 

9 11 

tive  conductivity  in  the  range  ~ 10  - 10  esu  [Kopecky  and  Kopecky, 

1971] • The  perpendicular  velocity  of  I50  m/s  required  by  (!+3)  for  this 
case  is  then  not  far  from  observed  downdraft  velocities  ~ km/s. 


The  question  of  whether  a steady  state  can  be  reached  is  also  relevant. 
We  shall  show  in  a subsequent  paper  that  for  a current  sheet  where  the 
thickness  and  perpendicular  scale  length  satisfy  -t  « r , the  magnetic 


resistive  decay  time  is  approximately 


2Jta^  ^r 


For  a radial  gas  inflow  the  gas  convection  time  is  clearly 


tconv  ~ r/v  • 
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Thus,  the  magnetic  Reynolds  number  for  a current  sheet  in  the  plane  of 
the  photosphere  is 


Rm  = tdiff/tconv  ~ 1 . (46) 

c 

We  note  that  the  self-consistency  relation  (43)  is  equivalent  to  R ~ 1 

m 

for  a current  sheet.  Even  for  values  -f,  r ~ 100  km  , Equation  (44) 
yields  decay  times  consistent  with  observation  [Stenflo,  1976].  The 
actual  lifetime  of  the  fluxules  is  probably  determined  by  MHD  considera- 
tions: e.g.,  how  long  can  the  low  pressure  core  be  maintained?  An  MHD 
instability  could  also  develop  to  destroy  the  dynamo  on  a time  scale 
related  to  (44). 


IV.  INHOMOGENEOUS  MAGNETIC  FIELD 


We  shall  restrict  our  attention  to  the  simplest  case  of  a vertical 

magnetic  field  B(x^)  which  is  the  same  in  both  slabs  and  also  treat 

2 = £ = E(x  ) . If  we  take  the  divergence  of  (12)  for  each  slab  and 

3D  1 

add,  we  obtain  using  (10) 

V • [se]  = - | V • [s  f e(Va  + Vb)]  . (47) 

The  field-aligned  current  density  is  thus 

•>«<*. *)  -£*  ■ [s  !«(vb  - ».>]  • <w) 

If  we  use  the  approximate  forms  of  (2)  and  (3)  for  the  dynamo  inequali- 
ties (4)  and  treat  n^  as  a constant,  Equation  (48)  can  be  cast  in 
the  form 

Jz(*,y)  vtjT,(VT«>  * vx[^i(1*a)(Vva>]|  < ^ 

this  result  and  (48)  reducing  to  (23)  for  V£  = VB  = 0 . For  the 
symmetric  case  we  have  treated,  the  polarization  field  need  not  be  com- 
puted to  obtain  J ; the  surface  currents,  of  course,  require  a solution 
for  E . Equation  (49)  gives  the  interesting  result  that  even  for  a 
constant  velocity  field,  a field-aligned  current  can  be  generated  through 

V X V ft  . , where  ft  . = eB/Mc  , when  a « 1 . 
ci  ci 
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V.  CONCLUDING  REMARKS 


I 

I 


We  have  considered  the  currents  generated  in  two  coupled  unipolar 
inductors.  For  simplicity  we  have  concentrated  on  the  case  of  a uniform 
magnetic  field  for  which  all  the  currents  may  be  computed  in  2 + 1 
Cartesian  geometry  with  relative  ease.  More  accurate  modelling  is 
needed  to  encompass  the  more  realistic  case  of  non-uniform  B . We  have 
indicated  in  Section  IV  how  gradients  in  magnetic  field  strength  can 
also  generate  field -aligned  currents  in  the  presence  of  neutral  wind 
flows . 

For  application  to  the  solar  atmosphere  we  have  considered  the  broad 
aspects  of  two  basic  categories  of  problems : 1)  in  the  presence  of  a 

strong,  large-scale  magnetic  field  (with  a subphotospheric  primary  current 
source)  can  photospheric  winds  distort  the  field  and  store  magnetic  energy 
in  the  corona,  to  serve  as  free  energy  for  a flare?  and  2)  can  the  photo- 
sphere generate  magnetic  fields  of  itself?  We  have  presented  some  evidence 
that  the  answer  to  both  questions  is  yes. 

In  regard  to  question  (l),  Heyvaerts  (19714-)  has  given  quantitative 
arguments  that  the  photosphere  can  be  a coronal  current  generator  using 
reasonable  values  for  parameters.  He  has  also  issued  the  caveat  that  for 
the  case  of  a thin  dynamo,  the  thickness  of  the  slab  is  important  in 
dimensional  analyses.  We  have  attempted  to  generalize  his  work  for  an 
arbitrary  gas  flow.  In  the  case  of  a bipolar  magnetic  arcade,  if  the 
dynamo  action  is  far  removed  from  the  polarity  inversion  line,  our  solution 
should  be  valid;  in  the  vicinity  of  the  inversion  line,  edge  effects  can 
occur  that  can  be  accounted  for  by  a modification  of  Equation  (29)  for  a 
semi-infinite  plane. 
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In  regard  to  question  (2),  we  believe  that  electrodynamically  there 
is  no  inconsistency  in  interpreting  small-scale  magnetic  flux  elements 
as  photospheric  dynamos.  The  magnetohydrodynamics,  of  course,  must  still 
be  done  to  get  an  overall,  self-consistent,  gas  flow  solution.  The  central 
problem  in  tying  everything  together  must  be  an  accurate  determination  of 
the  electrical  conductivity  tensor  a^(i=0=3)  for  each  solar  surface 
phenomenon  under  consideration. 

For  application  to  other  space  environments,  the  coupled  inductor 

model  can  provide  order  of  magnitude  estimates  for  many  effects  if  not 

just  valuable  insight.  The  important  point  is  that  the  model  attempts 

to  account  for  the  complete  current  circuit  and  in  doing  so  the  source- 

theoretic  approach  can  handle  a problem  where  local  dynamo  action  can 

produce  a global  effect  that  shows  up  a large  distance  away.  Field- 

theoretic  approaches  that  concentrate  locally  on  the  physics  in  a small 

volume  must  properly  take  account  of  boundary  conditions  far  away.  For 

instance,  the  influence  of  £.  on  a dynamo  in  slab  a is  manifest  through- 

fa 

out  (20) -(23).  With  appropriate  scale  factors  that  effect  the  transforma- 
tion from  a4^b  [viz..  Equation  (10)],  slab  b may  represent  a planet's 
ionosphere.  Slab  a can  be  used  to  describe  the  motion  of,  say,  the  iono- 
sphere of  Io  orbiting  Jupiter.  A generalized  Ohm's  law  may  be  constructed 
to  replace  (12)  for  slab  a and  the  model  may  be  able  to  provide  a crude 
description  of  the  complex  magnetosphere/ionosphere  interaction.  We  are 
pursuing  this  line  of  reasoning  in  the  hope  that  temporal,  electromagnetic 
induction  effects  may  be  incorporated  into  the  model. 
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DYNAMO  ACTION  IN  A THIN  SLAB 


by 

D.  D.  Barbosa 

Institute  for  Plasma  Research 
Stanford  University 
Stanford,  California  9^305 

FIGURE  CAPTION 

Figure  1.  Dynamo  action  in  two  thin,  coupled  slabs. 

(a)  Heyvaert's  (197*0  model  for  a bipolar  magnetic 
configuration.  Photospheric  convection  generates  a field- 
aligned  coronal  current  system  which  can  twist  or  shear 
the  field  and  provide  magnetic  free  energy  for  a flare. 

(b)  Two  resistive,  Cartesian  unipolar  inductors  connected 
in  series. 
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